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Abstract. 

Due to the presence of a large orbital angular momentum of the parton system 
produced at the early stage of non-central heavy-ion collisions, quarks and anti-quarks 
are shown to be polarized in the direction opposite to the reaction plane which is 
determined by the impact-parameter and the beam momentum. The global quark 
polarization via elastic scattering was first calculated in an effective static potential 
model, then using QCD at finite temperature with the hard-thermal-loop re-summed 
gluon propagator. The measurable consequences are discussed. Global hyperon 
polarization from the hadronization of polarized quarks are predicted independent 
of the hadronization scenarios. It has also been shown that the global polarization 
of quarks and anti-quarks leads also to spin alignment of vector mesons. Dedicated 
measurements at RHIC are underway and some of the preliminary results are obtained. 
In this presentation, the basic idea and main results of global quark polarization are 
presented. The direct consequences such as global hyperon polarization and spin 
alignment are summarized. 



1. Introduction 

We all know that spin is a basic degree of freedom of the elementary particles. Spin 
effects in high energy reactions usually provide us with useful information on the reaction 
mechanism and often give us surprises. Such effects have been studied intensively in 
high energy lepton-hadron, hadron-hadron, and hadron-nucleus collisions and lead to an 
active field of High Energy Spin Physics. In contrast, less study has been made in this 
direction in high energy heavy-ion collisions. One of the reasons might be that it would 
be very difficult or even impossible to polarize a heavy ion beam. Recent studies [1-3] 
have shown that hadrons can be polarized w.r.t. the reaction plane in high energy AA 
collisions with unpolarized beams. This is one of the places where we can study spin 
effects without a polarized beam and I hope that my talk can serve as an example which 
shows you, by looking at the spin degree of freedom, one can obtain some interesting 
information on the reaction mechanisms, even in heavy ion collision. 
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2. Global orbital angular momentum and shear flow 

We consider two colliding nuclei with the projectile of beam momentum pj„ moving in 
the direction of the z axis, as illustrated in Fig. [H The impact parameter h is taken 
as x-direction. The normal n;, oc pin x b of the reaction plane is taken as y. For a 
non-central AA collision, the dense matter system in the overlapped region will carry 
a global orbital angular momentum Ly in the direction —y. The magnitude of Ly is 
estimated using a hard spherical distribution for nucleus and is given in Fig. 2a. We 
see that —Ly is indeed huge and is of the order of 10^ at most 6's. 



1 X .V 




Figure 2. a (left), The global orbital angular momentum Ly of the overlapping system 
in a non-central AA collision at RHIC energy as a function of the impact parameter 6; 
b (right), The average orbital angular momentum of two neighboring partons separated 
by Ax = Ifm as a function of x/{Ra ~ fe/2) for different values of b/RA- 



Assuming that a partonic system is formed immediately after the initial collision, 
interactions among the produced partons will lead to the formation of a quark-gluon 
plasma (QGP) with both transverse (in x-y plane) and longitudinal collective motion. 
The existence of the global orbital angular momentum of the system discussed above 
implies a finite transverse (along ct) gradient of the longitudinal flow velocity. 

The initial collective longitudinal momentum can be calculated as the total 
momentum difference between participant projectile and target nucleons. Since the 
measured total multiplicity in AA collisions is proportional to the number of participant 
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nucleonsj4j, we can assume the same for the produced partons with a proportionahty 
constant c(s) at a given center of mass energy ^/s. Hence, the average collective 
longitudinal momentum per parton is given by, 



( h n _ dNP^Jdx - dNlJdx 

Pz{x,b; ^/s) is a monotonically increasing function of x. This can be seen more clearly 
by looking at the derivative dp^/dx, which is almost a constant for different x. From 
dpz/dx, we can estimate the average longitudinal momentum difference Apz between 
two neighboring partons separated by a transverse interval. On the average, the relative 
orbital angular momentum for two partons separated by Ax in the transverse direction 
is given by ly = —{AxYdpz/dx. E.g., for Au + Au at y/s = 200 GeV, c(s) ~ 45, ly for 
Ax = Ifm is shown in Fig. 2b. We see that ly is in general of the order of 1 and is larger 
than the spin of a quark. This implies that the effect can indeed be very significant. 

We emphasize that the results given in Fig. 2 correspond to the average over 
rapidity, and in principle, we should consider the distribution of such collective 
longitudinal momentum over rapidity. This kind of distribution has been discussed 
e.g. by Adil and Gyulassy in their study on jet tomography of twisted QGP [5], and 
can also be calculated e.g. using HIJING. If we take this distribution into account and 
consdier that partons with slightly different rapidities can also interact with each other, 
we will still obtain a non-vanishing vorticity for the interacting parton system. Here, 
the final relevant quantity should be the local (both in transverse separation x and in 
rapidity rj) derivative of the longitudinal momentum distribution. The average results 
presented above can only be serve as a guide for the magnitude of this effect. In the 
following, we will discuss the polarization effects caused by such local vorticity. 



3. Global quark polarization w.r.t. the reaction plane 

To see whether the local orbital angular momentum between the neighboring partons 
in a QGP can be converted to quark polarization via parton scattering, we consider 
quark scattering at fixed impact parameter. For definiteness, we consider a non- 
identical quark-quark scattering gi(Pi, Ai) + q2{P2, ^2) li^Pzi A3) + q2{Pi, A4), where 
Pi = {Ei,pi) and Aj denote the 4-momentum and spin of the quark respectively. We 
start with the usual cross section in momentum space, 

where M. {Q) is the scattering amplitude in momentum space, Q = P^ — Pi = P2 ~ -P4 is 
the 4-momentum transfer, Cgq = 2/9 and F are the color and flux factors, respectively. 
The cross section in impact parameter space is obtained by making a two dimensional 
Fourier transformation of the transferred transverse momentum qr, i.e., 

d^XT IQF.^f^J (27r)2(2vr)2 A(g^) a*(^^) ' ^ ^ 
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where A^(^t) and A^(fcT) are the scattering matrix elements in momentum space with 4- 
momentum transfer Q = {Q,q) and K = (0, k) respectively, A(g^) = ^ {Ei + £'2)!^ + Qzl 
is a kinematic factor. The differential cross section can be divided into a spin- 
independent and a spin dependent part, i.e., 
dax3 da ^ dAa 

Parity conservation demands that they have the following form, 

da ^ , /—^ dAa ^ , ^ ^ \ a ^ / ^\ / \ 

F(a;T,Vs), -— = n • (xr x p )AF(a;r, Vs), (5) 



d'-XT ' ' d^XT 

where n is the polarization vector for in its rest frame. F(xt, VS) and AF(xt, V§) 
are both functions of xt = \xt\ and the energy \/§ of the quark-quark system. This 
is because, in an unpolarized reaction, the cross section should be independent of any 
transverse direction. For the spin-dependent part, the only scalar that we can construct 
from the vectors that wc have at hand is n ■ (px xt). 

We note that, Xt x p is nothing else but the relative orbital angular momentum 

— * 

of the gig2-system, i.e., Xt x P — I- We see from Eq.(5) that the cross section takes 

— * — * 

its maximum n is parallel to Z or — Z depending on whether AF is positive or negative. 

— * — * 

This corresponds to a polarization of quark in the direction I or —I. 

As discussed in last section, for A A coUisions with given reaction plane, the direction 

— * — * 

of the averaged / of the two scattered quarks is given. Since a given direction of / 
corresponds to a given direction of Xx, this implies that there should be a preferred 
direction of Xt over others at a given direction of b. The detailed distribution of Xt at 
given b depends on the collective longitudinal momentum distribution discussed above. 
For simplicity, we considered an uniform distribution of xt in all the possible directions 
in the half oxy-plane with x > 0. In this case, we need to integrate da/d'^xx and 
dAa/iPxT in the half plane to obtain the average cross section at a given 6, i.e., 

r+00 r+ao da r+co r+oo dAa 

a = dx dy — Aa = dx dy — -=^, (6) 
Jo J-00 d'^XT Jo J-00 d'^XT 

The polarization of the quark after one scattering is given by, Pq — Aa jo. 
3.1. Results under small angle approximation 

The calculations are in principle straight forward but in practice very much complicated. 
Hence, in Ref.[l], we have given an example by calculating them using a screened static 
potential model and in the "small angle approximation" . The results are given by, 

da 1 g^CT f cPqr d^kr e'^kT-qryxr 



d^XT^SPM 2 ; {2ny {27ry {q^ + ^iDiPx + fil) 



(7) 



dAa^ _ .g^cr f d^qr d'^kr {kr - qr) ■ (p x n)e^('=^-9T)-^^ 
SPM - ~ J (2^^ (2^ p'(g^ + /^|p)(4 + /^D) ■ 



^d'^x 

Carrying out the integrations over qx and kx, we obtain that, 

'^'^ ' '^OilcTj^KlillDXT), (9) 



d^XT^SPM ' (27r)^ 
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Kq{iidXt)Ki{^dXt)), 



(10) 



SPM " ■ 

where Jq and Kq are the Bessel and modified Bessel functions respectively. Carrying 
out the integrations in the half plane with x > and we obtained that[l], 

Pg = -nfiDp/2E{E + mg). (11) 

The result is very encouraging since it shows a quite significant negative polarization of 
the quark after one scattering. 

More accurate calculations should be made using QCD at finite temperature. The 
quark-quark scattering is described by a Hard-Thermal-Loop (HTL) re-summed gluon 
propagator |6], 



A^'^(Q) 



+ 



+ (a-l)- 



where Q is the gluon four momentum, a is a gauge fixing parameter, 
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(12) 

(13) 

(14) 
(15) 



where uj = Q-U, Q = Q- luU, = -Q'^, g^y = g^y - U^Uu, x = uj/q, 
= g'^{N(. + Nf /2)T'^ /3 is the Debye screening mass, U and T are respectively the 
fiuid velocity and temperature of heat bath. In this framework, we have, 

Miqr) = ux,{Pi + Q)i,uxAPi)^'"'{Q)uxAP2 - QhuUx,{P2), (16) 
We work in the center of mass frame, and in Feynman gauge, so we have. 



A^'^(g) 
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+ 



(17) 



q^ q^ + fin^ 

The HTL gluon propagator needs to be regularized. We do this by introducing a non- 
perturbative magnetic mass fim ~ 0.255^jNj2g'^T [7] into the transverse self-energy. 
For simplicity, we consider only the longitudinal momentum distribution of the partons 
in QGP and have, in cm. frame, = (1, 0, 0, 0). 

Under the "small angle approximation", we obtain that. 



da 

d^XT 



9^Cqq 



dAa 

d^XT 



d^qT d'^k-T 
(27r)2 (27r)2' 

■9%q 
16 

X { 
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(19) 



We see that the difference between the results obtained with HTL propagator and those 
in the static potential model is the additional contributions from the magnetic part. 
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3.2. Numerical results 

The expressions for the cross sections without "small angle approximation" are quite 
complicated. I will not present them here. Interested readers are referred to [3]. We 
have carried out the integrations numerically and obtained the preliminary results as 
shown in Fig. 3a. For comparison, we show the results together with those obtained 
under the small angle approximation in Fig. 3b. 
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Figure 3. a (left), Preliminary results for the quark polarization -Pq as a function of 
\/l/T obtained using HTL gluon propagator, b (right). Comparison with the results 
obtained under small angle approximation (dashed line) and that obtained using static 
potential model under small angle approximation (dotted line). 



The results in Fig. 3 show that the polarization is quite different for different V^/T. 
It is very small both in the high energy and low energy limit. However, it can be as 
high as 20% at moderate v^/T. At RHIC, we can only give a rough estimation of the 
ratio \fI/T ~ Ap^/T which should be between 0.1 and 2. We see that in this range the 
polarization can be quite significant but can also be only of a few percent. 

From Fig. 3b, we see also that, at RHIC energy, small angle approximation is not 
a good approximation pp. We have to rely on the numerical results obtained using HTL 
gluon propagator without small angle approximation. More detailed study including 
transverse flow of the partons are also underway. 

4. The measurable consequences 

The global polarization of quarks and anti-quarks in QGP should have many measurable 
consequences for the hadrons after hadronization. The most direct ones are the 
polarizations of the spin non-zero hadrons. Data[8-10] from LEP on e"^e~ ^ Z ^ h + X 
tell us that polarization of the quark or anti-quark can indeed be transferred to final 
hadrons via hadronization. We now present the results for hyperons and vector mesons 
from the global quark polarization discussed above in the following. 
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^.i. Hyperon polarization 

For hyperons produced via recombination qqq —* H, we obtain, 

Pa = Ps-, = (4P, -Ps- 3PsP;)/{3 - AP.Ps + P^); (20) 

Ph = (4P. - P, - 3P,P,2)/(3 - 4P,P, + P2); (21) 

We see in particular that Ph = Pq for all the = A, S and S if P^ = Pq. 

For those produced via the fragmentation g ^ + X, we compare with the 
longitudinal polarization of hyperons in e^e~ ^ qq ^ + X which has been 

measured [8] and can be explained [11] by assuming that polarized hyperons contain the 
initial polarized leading quark in its SU(6) wave-function. Similar calculations lead to, 

Pa = nsPs/{ns + 2/,), P^ = (4/,P, - n,P,)/3(2/, + n,), (22) 
Ph = (4n,P, - /,P,)/3(2n, + /,), (23) 

where and fs are the relative s-quark abundances to u and d in QGP and 
fragmentation. We see in particular that Ph = Pq/'i if Ps = Pq and fs = Ug. 

In dependent of the hadronization mechanisms, we expect: (1) Hyperons and their 
anti-particles are similarly polarized; (2) Different hyperons are also similarly polarized. 
(3) The polarization vanishes in central collisions and increases with b in semi-central 
collisions. (4) It should have a finite value for small pt and central rapidity but increase 
with rapidity and eventually decreases and vanishes at large rapidities. 



4-2. Vector meson spin alignment 

The polarization of a vector meson V is described by the spin density matrix 
where the diagonal elements pY^, and p^i^i are the relative intensities for the spin 
component m of to take 1, 0, and —1 respectively. Pqq can be determined by measuring 
the angular distributions of the decay products. Furthermore, unlike the polarization of 
hyperons, pj^ does not know the direction of the reaction plane. Therefore, one cannot 
measure the sign of the quark polarization through . On the other hand, one does 
not need to determine the direction of the reaction plane to measure p^. 

For vector mesons produced in quark recombination mechanism, we obtain, 

P'ot^ = (1 - n')/(3 + P;), PoT^^^^^ = (1 - n^.)/(3 + PqPs). (24) 

We see in particular that p^ < 1/3 if is produced in this hadronization scenario. 

For the fragmentation of a polarized quark q^ —>■ V + X, we again compare with 
e"'"e~ ^ ^ qq ^ V + X , and obtain, 

P (frag) 

"3-/5P2' "n. + /.3-/3P2 + n. + /,3-/3Pr ^''^ 

The parameter P ^ 0.5 was obtained [12] by fitting the e~^e~ data[9,10]. We see that, 
in this case, poo > 1/3. We also see that, in both hadronization scenarios, vector meson 
spin alignment is a P^ effect. This has the advantage discussed at the beginning of this 
sub-section but also the shortage that it may be very small. 
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We are happy to know that dedicated efforts have been made in measuring such 
effects at RHIC. See e.g. [13,14]. 

5. Summary and outlook 

In summary, we have shown that produced partons have large local relative orbital 
angular momentum in non-central AA collisions at high energies. Parton scattering 
with given relative orbital angular momentum can polarize quarks along the same 
direction due to spin-orbital interaction in QCD. Such global quark polarization has 
many measurable consequences and the measurements on such effects might open a new 
window to study the properties of QGP in high energy AA collisions. 
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